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Using the bosonization and level spectroscopy methods, we study the ground-state phase di- 
agram of a XXZ antiferromagnet on a railroad-trestle lattice with asymmetric leg interactions. 
It is shown that the asymmetry does not change the dimer/Neel transition line significantly, 
which agrees with the expectation based on a naive bosonization procedure, but it does change 
the dimer/spin- fluid transition line. To understand this observation, we analyze eigenvectors of 
the ground state, dimer excitation, doublet excitation and Neel excitation, and find that only 
the doublet excitation is affected by the asymmetric interaction. 

KEYWORDS: railroad-trestle lattice, sawtooth lattice, XXZ model, exact singlet dimer state, level spec- 
troscopy method, bosonization 
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1. Introduction 

Spin-1/2 Heisenberg antiferromagnets on a railroad- 
trestle lattice have been investigated intensively, be- 
cause they show interesting quantum phase transitions 
and critical phenomena. 1-5 Haldane 1 and Nakano and 
Fukuyama 2,3 discussed the phase diagram by using the 
bosonization and renormalization group methods. Tone- 
gawa and Harada numerically estimated the dimer-Neel 
transition point. 4 It is, however, known that a quan- 
titative estimation of the critical point is difficult be- 
cause of the nature of the Berezinskii-Kosterlitz-Thouless 
(BKT) transition. To overcome this difficulty, Nomura 
and coworkers developed the level spectroscopy (LS) 
method, in which the transition point between the dimer 
and Neel (spin- fluid (SF)) phases is estimated by the in- 
tersection of energies of the dimer and Neel (doublet) 
excitations, and they succeeded in obtaining a quantita- 
tive result of the phase diagram for the XXZ model on 
a railroad-trestle lattice. 5-9 

The sawtooth-lattice Heisenberg antiferromagnet has 
been studied by several authors because of the following 
two reasons: (i) this system is one-dimensional counter- 
part of the Kagome lattice and shows a two-stage entropy 
release as a function of temperature, 10 ' 11 and (ii) the 
sawtooth lattice is realized in YCUO2.5. 12 In this paper, 
we consider a model that contains both the sawtooth 
and railroad-trestle models as special cases, and study 
the ground-state phase diagram in a unified way. 

We introduce the model Hamiltonian 



H = Y,Hi + i + aJ2 {1/2 + {-ITS} h iti+2 , 



(1) 



with h h3 = SfS* + SfSV + AS?S*. (See Fig. 1(a).) Here, 



L denotes the total number of spins, and the periodic 
boundary condition Sl+i = S± is assumed. We also con- 
fine ourselves to < a < 1, — 1 < A and < 8 < 1/2. 
This model is reduced to the railroad-trestle model and 
sawtooth model, respectively, at 6 = and 1/2. 13 The 
quantitative ground-state phase diagram of the railroad- 
trestle model with A > was obtained in ref. 6 for 



the first time, and the extension to the A < region 
was made in ref. 9, where a direct transition between 
the dimer state and the ferromagnetic state was found. 
Here, we study how the asymmetry parameter 8 affects 
the phase boundaries. 




q(1/2-«) 



1 + 2 



Fig. 1. Graphical representations of (a) the asymmetric railroad- 
trestle model and (b) a four spin problem (see text). In (b), each 
dashed line represents ^^+2 or — and the thin solid line 

represents a reflection plane. The hatched spin pairs are used to 
define the two-spin basis set. 



This paper is organized as follows: in §2, we derive 
some basic properties of the present model. We show 
that the model has an exact dimer singlet ground state 
at a = 1, and a naive bosonization procedure predicts 
that critical lines in the ground-state phase diagram do 
not depend on 8. In §3, we describe our procedure of nu- 
merical calculations based on the LS method. We present 
the ground-state phase diagram in the a — A plane in §4, 
and discuss the ^-dependence of the transition lines in 
§5. We summarize our results in §6. 

2. Basic Properties of the Model 

At a = 1, the system has an exact dimer singlet ground 
state, independent of 8. To see this, we first write the 
singlet dimer state as 



1 



and define 



[i,j]=j={\Ui)-\Uh)}, 



$i(L) = [1, 2] [3, 4] •••[£-!, L], 



(2) 
(3a) 
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$ 2 (L) = [2, 3] [4, 5] •••[£, 1]. (3b) 
We rewrite eq. (1) with a = 1 as 

H a=1 = H a=liS=0 + S ^(-1)^^+2- (4) 

i 

The first term, if a =i,5=o, is simply the railroad-trestle 
model at the Majumdar- Ghosh point, and thus, it has 
doubly degenerate exact dimer singlet ground states, 
<£>i(L) and <3?2(£). 14 ~ 16 In order to show that the dimer 
singlet state is an eigenstate of eq. (4), it is convenient 
to use the dimer basis set and the reflection operation 
about the thin solid line shown in Fig. 1(b). The two- 
dimer singlet state [i, i + + 2, i + 3] has even parity 
under both the reflection and spin inversion operations. 
The operation of hi^ +2 — ^i+i,i+3 to a two-dimer state 
changes the parity under the reflection operation. Thus, 
(hij+2 - h i+ ij +3 )[i,i + l][i + 2, i + 3] is an odd-parity 
state of the reflection. Among all two-dimer states with 
S*t ot = 0, there are two odd-parity states of the reflec- 
tion, [i,i+ l](i + 2,i + 3) and l)[i + 2, i + 3] with 
{hi) = {\hij) + |.UTj}}/v / 2- However, these two states 
have odd parity under the spin inversion. Thus, we ob- 
tain {h l:i+ 2 — hi+i,i+3)[i,i + + 2,i + 3] = 0, which 
shows that $i(£) and $ 2 {L) are eigenstates of H a= \, 
even with 8 ^ 0. The energy eigenvalue of these singlet 
dimer states is given by 

A + 2 T 



are 



Edh 



(5) 



In order to prove that the dimer singlet state is a ground 
state, we rewrite the Hamiltonian in eq. (4) as the sum 
of spin plaquette parts: 



L/2 



(6) 



1=1 



with 



hf aq = h 2l , 2l+1 + (1/2 - 5){h 2 i-i,2i + h 2 i-i,2i+i) 

+ (1/2 + 6)(h 2 l+l,2l+2 + &2J.2J+2)- (7) 

The minimum eigenvalue of hf laq is -(2 + A)/4 for A > 
— 1/2 and 3A/4 otherwise. Therefore, the ground-state 
energy E g of H a= i satisfies 



Eg> 



4l 



for A > 
for A < 



-1/2 
-1/2 



(8) 



The lower bound of E g for A > —1/2 agrees with -Edimen 
which proves that the dimer singlet state is a ground 
state. For A < —1/2, the ferromagnetic state, whose 
energy is -Ef or ro = is a ground state. 

We turn to the general case with a / 1. To obtain 
insight into the ground-state phase diagram, we bosonize 
the Hamiltonian in eq. (1), which leads to 



cos - 



(9) 



H=lj dx[ A m f + B { ^f]-cf {2nar 

where a is a lattice constant. The coefficients A, B, and 
C are, in general, functions of A, a and 8. In the vicinity 
of A = and a = 0, the explicit forms of A, B, and C 



B=1 _ A_a(2-A) 



C =a A 



7T 2lT 

a(2 + A) 



(10a) 
(10b) 
(10c) 



Note that A, B, and C do not depend on 5, which is 
because 5 Eleven and S^2 ieodd h iti+2 in eq. (1) 

cancel each other out in the long-wavelength limit. This 
fact suggests that the phase boundaries in the a — A 
plane do not change when 5 varies from the railroad- 
trestle point 8 = to the sawtooth point 8 = 1/2. 

At the Hciscnberg point, A = 1, the phase boundary 
between the dimer and Neel phases has already been de- 
termined for both the railroad-trestle and sawtooth mod- 
els by using the LS method. Previous studies estimated 

a dimcr/Nccl (A = ^ = Q) = Q ^ fof ^ raUroad _ 

trestle model 5 and a dimcr/Ncci (A = 16 = 1/2) = q 4374 

for the sawtooth model, 17 which indicates that the re- 
lation a dimcr/Nccl (A = ^ = Q) = ^dimcr/Nccl ( A = 

1,8 = 1/2) suggested by the bosonization method ap- 
proximately holds. 

In this paper, we calculate the ground-state phase di- 
agram of the asymmetric railroad-trestle model by using 
the LS method and study how the phase boundaries in 
the a — A plane depend on 8. As a result, we find that 
the 5-dependence of the dimer-Neel critical line is small, 
but that of the dimer-SF critical line is large. 

3. LS Method 

We now describe our calculations based on the LS 
method. 5 " 8 In Fig. 2, we show the excitation energies 
of Neel, dimer, and doublet excitations as functions of 
a in a finite-size cluster with L = 20. The intersection, 
a c (L; A = 0.5, 8 = 0), between the dimer and doublet ex- 
citation energies in Fig. 2(a) is interpreted as a dimer-SF 
transition point in the L = 20 cluster. Also, the intersec- 
tion, a c (L; A = 2, 8 = 1/2), between the dimer and Neel 
excitation energies in Fig. 2(a) is interpreted as a dimer- 
Neel transition point. The size dependence of a c (L; A, 8) 
is expected to be 5,6 

a c (L; A, 8) = a c (oo; A, 8) + Const, x L~ 2 . (11) 

Our calculated data of a c (L; A, 8) with L = 12, 16, • • • , 
28 are plotted against L~ 2 in the insets in Fig. 2. We 
can confirm from this figure that our data follow the 
size dependence in eq. (11) and that the extrapolation 
to L = 00 is satisfactory. 

Next, we calculate the central charge c for verifying our 
calculations. According to conformal field theory, we can 
estimate c using the following expression for the ground- 
state energy: 



E a {L) = e L- 



TTVC 

~6L ' 



(12) 



with the spin-wave velocity v = LAE(k = 2tt/ L)/2ir. 
Figure 3(a) shows the (5-depcndcncc of c for (8, A) — 
(0,0) and (0.5,2,5). As expected, it is found that c = 1 
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in the SF phase {a < 0.6474) for (6, A) = (0, 0) and at 
the dimer-Neel transition point (a = 0.7048) for (S, A) = 
(0.5,2.5). We also calculate c along the dimer-SF and 
dimer-Neel critical lines, where c = 1 is expected. The 
result is shown in Fig. 3(b). Our result is consistent with 
c = 1 , although a small error is observed. 
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point 8 = 1/2. The pronounced <5-dependence of the 
dimer-SF critical line contradicts the naive expectation 
based on the bosonization. 




Fig. 4. The ground-state phase diagram in the a — A plane for 
<5 = 0, 1/4 and 1/2. 
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Fig. 2. Excitation energies of Neel, dimer, and doublet excita- 
tions as functions of a for (a) 5 = 0, A = 0.5 and (b) 5 = 1/2, 
A = 2, where the system size is L = 20. 





+ 8-0 ;A=0 
X 5=1/2;A=2. 

- 8=0 trans 

8=1/2 trans 


tion point 
ition point 




0.0 0.2 0.4 0.6 0.8 1.0 



(b) 1.2 









+ 8=0 






X 8=1/2 




- x J ¥ ¥ * * ; 


* ¥ If 


dimcr/SF 


dimcr/Nccl - 







-0.5 0.0 0.5 1.0 

A 



0.5 0.0 

1/A 



Fig. 3. (a) Plot of the central charge c against a for fixed values 
of 5 and A. (b) Central charge c on the critical lines as a function 
of A. 



To see the (5-dependence of the transition lines in more 
detail, we define 



a'(A,5) = a c (A,<5)-a c (A,0) 



(13) 



and show a'(A, 1/4) and a'(A, 1/2) in Fig. 5. We find 
that a' (A, 1/2) for the dimer-Neel and ferro-dimer tran- 
sitions is at most <~ 10~ 2 . For the dimer-SF transition, 
a'(A, 1/2) increases as A decreases, and it amounts to 
<~ 10 _1 at approximately A = —0.5. Roughly speaking, 
the (5-dependence of the dimer-SF transition is about ten 
times as large as that of the dimer-Neel and ferro-SF 
transitions. 




4. Ground-State Phase Diagram 

Putting together all our calculated data, we obtain the 
phase diagram in the a — A plane shown in Fig. 4, where 
the phase boundaries for 6 = 0, 1/4 and 1/2 are pre- 
sented. Our result for 6 = reproduces the ground-state 
phase diagram of the railroad-trestle model previously 
obtained by Nomura and Okamoto 6 and Hirata and No- 
mura. 9 We find in Fig. 4 that the dimer-Neel critical line 
and the first-order phase transition line to the ferromag- 
netic state are almost independent of 6. On the other 
hand, the dimer-SF critical line depends on S, and this 
dependence is pronounced at approximately A = —0.5. 
As a result, the region of direct transition between the 
dimer and ferromagnetic states 9 becomes narrower when 
5 increases, and this region vanishes at the sawtooth 



Fig. 5. Plot of a' (A, 8) = a c (A, 8) - a c (A, 0) as a function of A 
for 5 = 1/2 and 1/4. The markers are the results in the thermo- 
dynamic limit and the lines are those for L = 4. 



5. Discussion 

In this section, we study the eigenvectors for the 
ground state, dimer excitation, Neel excitation, and dou- 
blet excitation to discuss why the (5-dependence of the 
dimer-SF transition line is much larger than those of the 
dimer-Neel and ferro-SF transition lines. 

We begin by representing H graphically and defining 
a reflection operation I and a translation operation t, 
as shown in Fig. 6. Note that I and t do not change H 
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only when 6 = 0, but change H when 6 0. We also 
define the spin inversion operation as T. Hereafter, we 
denote the eigenvalues of T, i, and I as T, t = e lk , and 
/, respectively. 




Fig. 6. Another graphical representation of the Hamiltonian in 
eq. (1). We denote the translation operation as t and the reflec- 
tion operation about the horizontal thin solid line as I. 



Table II. Symmetrized basis functions with total S^ ot = 1,2 for 
the L = 4 cluster. 



(S t z ot ,k,I) 


basis function 


(1,0,1) 


(|*0*+) + \t+t ))/V2 


(1,T,-1) 


(\st + ) + \t+s))/V2 




(\st+) - \t+s) - i\t t+) + i|*+* >)/2 


(L-f.*) 


(\st+) - \t+s) + i\t t+) - i|*+* >)/2 


(2,0, 1) 


lM+> 



Table III. Total numbers of basis functions in the subspaces rel- 
evant to the determination of transition lines. 



state 


(T orS* ot ,fc,/) 


L = 4 


8 


12 


16 


ground 


(1,0,1) 


2 


7 


35 


257 


state 


(l,ir,-l) 





1 


13 


175 


dimer 


(l,ir,l) 


1 


5 


29 


239 


excitation 


(1,0,-1) 








9 


158 


Neel 


(-1,7T,-1) 


1 


4 


27 


230 


excitation 


(-1,0,1) 





1 


15 


183 


doublet 




1 


5 


38 


375 


excitation 


(1,0,1) 


1 


5 


38 


375 



When 5 = 0, it is known that the ground state is 
in a S£ ot = subspace with (T,k,I) = (1,0,1), the 
dimer excitation is in the subspace with (l,7r, 1), and 
the Neel excitation is in the subspace with (— l,7r, — l). 6 
Operation of ^(("l)'^^ to basis functions with 
(T,k,I) = (1,0,1), (l,7T,'l) or (-1,tt,-1) yields ba- 
sis functions with (T,k,I) = (l,7r, — 1), (1,0,-1) or 
(— 1, 0, 1), respectively. Thus, for 6 > 0, the ground-state 
wave function, \^ gs ), the wave function of the dimer exci- 
tation, |vE'dimor), and that of the Neel excitation, I^Ncci), 
are written as 

|$ g s) H0i,o,i) cos 0gs + l0i,7r,-i)sin0 gs , (14a) 

|$dimer) = |01,7r,l) COS $dimcr + |01,O,-1> Sm6» dimcr , (14b) 
|*Ncol) = |0— 1,7T, — 1> COS N <Sel + |^-l,0,l) Sin 6»Neel,(14c) 

where \4>T,k,i) represents unit vector in the (T, k, I) sub- 



space of |$ gs ), |$d 



or |$Ncei), respectively, and we 



choose < 9 < tt/2. (We hereafter call 9 the "mixing 
parameter".) For the doublet excitation, it is known 
that the symmetry of the doublet excitation is charac- 
terized by (5 t z ot , k, I) = (1, 7r, — 1) when 6 = 0. Operation 
of ^i(~^)'' i i,'+2 to those basis functions yields basis 
functions with (1,0,1). We can write the wave function 
of the doublet excitation for 6 > as 

|$doublct) = |<Pl,7T,-l) COS^doublct + |Vl,0,l) Sin ^doublet, 

(15) 



Table I. Symmetrized basis functions with S5 



: for the L ■■ 



cluster. The eigenvalues of T, i and I are denoted as T, e lk and 
respectively. 



(T,k,I) 


basis function 


(1,0, 1) 


(\ss) + \t t ))/V2 

(\ss) - |*o*o) - \t+t-) - \t-t+))/2 


(1,T,1) 


(|ss> - |* *o> + |*+*-> + l*-*+>)/2 


(-l,w,-l) 


(|st > + |tos»/v^ 


(-1.5.*) 


(|st > - \t s) - i\t+t-) + i\t-t+))/2 


(-1.-5.*) 


(\st ) - |t s> + i\t+t-) - i\t-t+))/2 



where \tps* ,k,i) is unit vector in the (5 t z ot , k, I) subspace 

Of |$doublet)- 

Here we consider the case of L = 4. It is known that 
LS analysis gives good results even for small systems; 
thus, we expect that the study of the L = 4 cluster is a 
good starting point for understanding the ^-dependence 
of the transition lines. The symmetrized basis functions 
are tabulated in Tables I and II, where the basis functions 
are expressed by dimcr states, \s) = -^(| tl) — I IT)), 

l*o) = 775G Tl) + I IT)), \t+) = I TT), and |t_) = I || 
). For example, a four-spin state \sto) means that two 
spins on the first and second sites are in the singlet dimer 
state and those on the third and fourth sites are in the 
triplet dimer state with S^ ot = 0. Note that there are no 
basis functions with (T,k,I) = (1, n, — 1), (1, 0, — 1), or 
(—1,0, 1). (See also Table III, where the total numbers 
of basis functions in each of the subspaces are listed.) 
Therefore, we obtain 

#gs = #dimer = #Nccl = 0, (16) 

which means that the energy eigenvalues of the ground 
state, dimer excitation, and Neel excitation do not de- 

n c- i dimcr/Nccl i SF/fcrro 

pend on o, and thus a c and a c are com- 

pletely independent of S in the L = 4 cluster. For the 
doublet excitation, however, there is a basis function in 
the (S* ot = 1, k = 0,1 = 1) subspace. Therefore, ^doublet 
does not vanish. Using the basis functions in Table II, we 
obtain 



^doublet = arctan 



v/l + M(l-A)P 



(17) 



a5(l - A) J 

which means that Q * mcr / SF depends on 6. An explicit 
calculation gives 

2 



a 



dimcr/SF 



■ (18) 

y/(3 + A) 2 -45 2 (A- l) 2 
Using the results obtained above, a'(A,6), which is dc- 
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fined in eq. (13), for the L = 4 cluster is calculated as 



2 

3+A 



L 3 + A \ 





- 1 



for dimer/SF 



otherwise 

(19) 

The present analysis of the L = 4 cluster shows that 
the pronounced 5-dependence of a' (A, 5) for only the 
dimer-SF transition is because only the mixing parame- 
ter doublet for the doublet excitation can become finite. 

In Fig. 5, we show a' L=4 (A,6) for S = 1/4 (dashed 
line) and 1/2 (solid line). We find that a' L=4 (A,S) 
shows a similar characteristic behavior to a' L=OQ (A, S). In 
particular, a' L=4 (A,S) agrees well with a' L=oc (A, 5) for 
—0.5 < A <^ 0.5. However, a' L=0O (A 7 5) for Neel-dimer 
and ferro-SF transitions is small but finite, although 
a' L _ 4 (A, S) = 0. As shown in Table III, basis functions 
with {T,k,I) = (1,tt,-1), (1,0,-1), and (-1,0,1) ap- 
pear for larger clusters, so it is possible to be 9 gs , #dimer, 
#N6ci 7^ 0. Even in such larger clusters, we expect that 
these values remain smaller than that of the doublet ex- 
citation. We confirm this fact numerically for a larger 
cluster (L = 12) below. 

We now study the mixing parameters for the L = 
12 cluster. Here, (A, a) is fixed to be (2.0,0.7) and 
(—0.5,0.9), where the former and latter are chosen to 
be near the dimer-Neel and dimer-SF transition lines, 
respectively. We calculate the mixing parameters as func- 
tions of 6. The result for (A, a) = (2.0,0.7) is shown in 
Fig. 7(a). We find that all the mixing parameters are fi- 
nite for S > and these values increase as S increases. 
Although #Neci is comparable to ^doublet, these values are 
very small, 0.04 x ir/2 at the most. The smallness of the 
mixing parameters results in the dimer-Neel transition 
line being almost independent of 5. Next, we show the 
result for (A, a) = (-0.5,0.9) in Fig. 7(b). As expected, 
we find that ^doublet is much larger than the other val- 
ues, and it amounts to 0.17 x 7r/2 at 5 = 0.5, which leads 
to the pronounced (^-dependence of the dimer-SF transi- 
tion line. Thus, we conclude that the feature extracted 
from the L = 4 cluster analysis, the 5-dependence of the 
dimer-SF transition line due to the mixing of the doublet 
excitation, also holds for larger clusters. 

6. Summary 

The ground-state phase diagram of the XX Z model on 
the railroad-trestle lattice with asymmetric leg interac- 
tions has been investigated. We have proven the existence 
of an exact singlet dimer ground state at a = 1, and 
have studied how the phase boundaries depend on the 
asymmetry parameter 5. The naive bosonization method 
predicted that the asymmetry does not affect the critical 
lines, but this prediction has turned out to be an oversim- 
plification. By using the LS method, we found that the 
dimer-Neel critical line and the ferro-SF first-order tran- 
sition line are almost independent of S, and the dimer-SF 
critical line depends on S, which is because the asymmet- 
ric interaction affects the doublet state, but it does not 
affect the ground state, dimer excitation or Neel exci- 
tation significantly. This may be interpreted as follows: 
the doublet excitation is a type of spin- wave state in the 



(a) A=2.0,<x=0.7,L=12 
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Fig. 7. Mixing parameter against 5 for (a) A = 2.0, a = 0.7 
and (b) A = —0.5, a = 0.9. The system size is L = 12. 



quasi long-range ordered state, and a spin wave, which 
is a local defect propagating in the quasi long-range or- 
dered state, can detect the short-wavelength part of the 
Hamiltonian, such as the asymmetric interaction. 

After completion of this work, we became aware of 
studies by Sarkar and Sen 18 and by Capriotti et al. 19 con- 
cerned with Heisenberg antiferromagnets on the asym- 
metric railroad-trestle lattice. Sarkar and Sen produced 
an operator product expansion of the terms in the 
bosonized Hamiltonian from the asymmetric part in 
leg interactions, and showed that these terms cancel 
each other. Capriotti et al. numerically calculated the 
dimer/SF transition point as a function of 5 at the 
Heisenberg point, and they found that the ^-dependence 
is very small. 
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